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A brief overview of our current understanding of abundance and properties of dark energy 
and dark matter is presented. A more focused discussion of supersymmetric dark matter 
follows. Included is a frequentist approach to the supersymmetric parameter space and 
consequences for the direct detection of dark matter. 

1 The Energy Density Content of the Universe 

The overall composition of the Universe can be conveniently described by the density parameter, 
Q,, defined as the average energy density of the Universe, p, relative to the critical density needed 
for a spatially flat Universe, pc- One of the Einstein field equations leads to the expression for 
the expansion rate of the Universe, which we characterize by the Hubble parameter, 

^ = [r) "i?^+3' 

where R{t) is the cosmological scale factor and k is the three-space curvature constant {k = 
0, +1, —1 for a spatially flat, closed or open universe). A is the cosmological constant which is 
assumed here to contain all contributions from the vacuum energy density. One can define a 
critical energy density pc such that p = pc for fc = 

Pa = SHySnGN- (2) 

In terms of the present value of the Hubble parameter this is, 

Pc = 1.88 X 10"2^/lo^gcm-^ (3) 

where 

ho = i/o/(100kmMpc"^s-i). (4) 
The cosmological density parameter is then defined by 

n = -^. (5) 

Pc 

The composition of the Universe can be expressed by breaking down the density parameter 
into separate contributions, 

n = nr + n,n + ^A, (6) 
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for contributions from radiation, matter and a cosmological constant /vacuum with JIa == 
A/3H^. The contribution to fir from the cosmic microwave background (CMB) is small, of 
order 10^'*. Precise determinations of the matter and vacuum contributions to are obtained 
from the detailed power spectrum of CMB anisotropics as measured by WMAP [T]. When 
combined with other measurements such as high redshift supernova type la data [2] and baryon 
acoustic oscillations [3], one finds 

/lo = 0.71 ±0.01 17o = 1.006 ±0.006. (7) 

WMAP data alone is sufficient for determining the individual contributions to ft of 

rJ™/i^ 0.133 ±0.006 r^A = 0.74 ±0.03. (8) 

The matter content of the Universe can be further broken down as WMAP also determines 
the baryon density of Universe [T] 

r^s/ig = 0.0227 ±0.0006. (9) 

The contribution to f2 in neutrinos lies in the range 

0.0005 < n^hl < 0.0076, (10) 

where the lower bound is obtained from the requirement of finite neutrino masses from oscilla- 
tion data and the upper bound is again derived from WMAP data in conjunction with other 
large scale structure data. 



2 Dark Energy 

The biggest surprise of all of the recent determinations of contributions to Q must be the 
realization that there is a substantial contribution from dark energy, namely that JIa 0. The 
WMAP value for Q\ is moreover consistent with determinations from supernovae data and 
baryon acoustic oscillations. When all data are used, one finds I^a = 0.726 ± 0.015. 

But now a bigger question arises: What is the physical nature of the dark energy? Different 
possibilities can be distinguished by their equation of state characterized hy w = p/p. The 
equation of state parameter for radiation is simply Wr = 1/3, whereas for matter, it is Wm = 0. 

The simplest solution for the dark energy remains either a cosmological constant or a con- 
stant vacuum contribution to the energy density with an equation of state, w = —1. This is 
indeed consistent with the central value determined by WMAP, which finds [1] 

-0.33 < 1 + wo < 0.21, (11) 

for the value of w today (at 95 % CL). The numerical value for A, however, is extremely small, 
and when written as a dimensionless constant (as GatA), it is of order lO"^^'^. This is the well 
known cosmological constant problem in cosmology [4] . 

There are, of course, other possibilities, the largest class of which is known as quintessence 
[5]. In this case, the dark energy may be a dynamical phenomenon described by an evolving 
scalar field. The energy density and pressure of a scalar field, ip, with potential, V{4>), can be 
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written as (neglecting spatial gradient terms) 



P=\^'' + V{cj,) (12) 

When the kinetic term is small compared to the potential, p k,V and p and we recover 

the constant solution with w = —1. In general, however, wq may differ from -1 and indeed may 
not even be a constant. Once again, WMAP (using supernovae and BAO data) place combined 
limits on w and its derivative with respect to redshift, w' , 

w = -1.06 ±0.14 w' = 0.36 ±0.62 (14) 

In short summary, we are left with the following puzzles regarding dark energy: 

• There is the question of fine-tuning associated with the cosmological constant problem. 
Namely, we expect several contributions to the vacuum energy density 

A = AguT + AeW + ^QCD ■ ■ ■ (15) 

where the various contributions listed arise from possible sources such as grand unified 
theories [GnAgut (10^^)^) , the Standard Model {GnAew ^ (lO'^'^Y), and QCD 
[GnAqcd - (10-20)4), yet sum to lO^i^s. 

• The coincidence problem. Here, we would like to understand why f2,„ is within a factor 
of a few of today. This is one of the issues addressed by quintessence models and may 
be probed in observations testing the possibility of variability in fundamental constants. 



3 Dark Matter 

From the quoted contributions to f2 in matter and baryons from WMAP, we can obtain the 
density of cold dark matter from the difference between the total matter density and the baryon 
density [1] 

f^cm/Zi^ = 0.1099 ±0.0062 (16) 

or a 2cr range of 0.0975 - 0.1223 for ilcDAih^- 

Evidence for dark matter in the universe is available from a wide range of observational data. 
In addition to the results from the CMB, there is the classic evidence from galactic rotation 
curves [6] , which indicate that nearly all spiral galaxies are embedded in a large galactic halo of 
dark matter leading to rather constant rotational velocities at large distances from the center 
of the galaxy (in contrast to the expected ^ 1/r behavior in the absence of dark matter). 
Other dramatic pieces of evidence can be found in combinations of X-ray observations and weak 
lensing showing the superposition of dark matter (from lensing) and ordinary matter from X-ray 
gas [7] and from the separation of baryonic and dark matter after the collision of two galaxies 
as seen in the Bullet cluster 0. For a more complete discussion see [5]. 

In addition to being stable (or at least very long lived), the dark matter should be both 
electrically and color neutral. Indeed, there are very strong constraints, forbidding the existence 
of stable or long lived particles which are not color and electrically neutral as these would become 
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bound with normal matter forming anomalously heavy isotopes. The limits on the abundances, 
relative to hydrogen, of nuclear isotopes [10], n/nn ^ 10^^^ to 10"^^ for 1 GcV < m < 1 
TeV. A strongly interacting stable relic is expected to have an abundance njuH ^ 10^ ^'^ with 
a higher abundance for charged particles. 

Unfortunately, there are no viable candidates for dark matter in the Standard Model. As 
baryons and neutrinos have been excluded, one is forced to go beyond the Standard Model, and 
here, I will focus on the possibilities which exist in the context of the minimal supersymmetric 
extension of the Standard Model (MSSM) [11]. In the MSSM, the lightest supersymmetric 
particle (LSP) is stable if R-parity (i? = — i3S+l+2s^ unbroken. There are several possibilities 
in the MSSM, specifically the sneutrino with spin zero, the neutralino with spin 1/2, and 
the gravitino with spin 3/2. However, a sneutrino LSP would have relatively large coherent 
interactions with heavy nuclei, and experiments searching directly for the scattering of massive 
dark matter particles on nuclei exclude a stable sneutrino weighing between a few GeV and 
several TeV [12] . The possible loophole of a very light sneutrino was excluded by measurements 
of the invisible Z-boson decay rate at LEP [13] . The gravitino is a viable candidate and often 
predicted in models based on supcrgravity |14l 115] . In this case, however, its probability for 
direct detection is negligible. 

There are four neutralinos, each of which is a linear combination of the i? = — 1 neutral 
fermions [TB]: the wino ^ the partner of the 3rd component of the SU{2)l gauge boson; the 
bino, B- and the two neutral Higgsinos, Hi and H2. The mass and composition of the LSP are 
determined by the gaugino masses. Mi and M2, the Higgs mixing mass term, /i, and the ratio 
of the two Higgs vacuum expectation values expressed as tan /3. In general, neutralinos can be 
expressed as a linear combination 

X = aB + pW^+-fHi+SH2. (17) 

The relic density of neutralinos depends on additional parameters in the MSSM beyond 
Ml, M2, n, and tan/3. These include the sfermion masses, and the Higgs pseudo-scalar 
mass, ruA- To determine the relic density it is necessary to obtain the general annihilation 
cross-section for neutralinos. In much of the parameter space of interest, the LSP is a bino and 
the annihilation proceeds mainly through sfermion exchange. 

In its generality, the MSSM has over 100 undetermined parameters. There are good argu- 
ments based on grand unification [17] and supcrgravity |18] which lead to a strong reduction 
in the number of parameters. I will assume several unification conditions placed on the super- 
symmetric parameters. In all models considered, the gaugino masses arc assumed to be unified 
at the GUT scale with value, mi/2, as are the trilinear couplings with value Aq. Also common 
to all models considered here is the unification of all soft scalar masses set equal to mo at the 
GUT scale. With this set of boundary conditions at the GUT scale, we can use the radiative 
electrowcak symmetry breaking conditions by specifying tan/3, and the mass, Mz, to predict 
the values of ^ and Higgs pscudoscalar mass, tua- The sign of fj, remains free. This class of 
models is often referred to as the constrained MSSM (CMSSM) [H] [20l [211 [22l [23] . In the 
CMSSM, the solutions for jj, generally lead to a lightest neutralino which is very nearly a pure 
B. 

I note that that while the name CMSSM is often used synonymously with mSUGRA, for 
minimal supcrgravity [181 124| . The latter however, has two additional constraints: 7TI3/2 ~ itiq 
and Bq = Aq — mg. The former sets the unification scale scalar masses equal to the gravitino 
mass. This condition often results in a gravitino LSP [14]. The latter condition sets a relation 
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between the supersymmetry breaking bilinear, trilinear and scalar mass terms. Because of this 
condition, tan/3 is no longer a free parameter, but must be solved for through the radiative 
electroweak symmetry breaking relations. 

In Fig. [l] an example of the renormalization group running of the mass parameters in the 
CMSSM is shown. Here, we have chosen mi/2 = 250 GcV, toq = 100 GeV, tan;3 = 3, = 0, 
and /i < 0. Indeed, it is rather amazing that from so few input parameters, all of the masses of 
the supersymnictric particles can be determined. The characteristic features that one sees in 
the figure, are for example, that the colored sparticles are typically the heaviest in the spectrum. 
This is due to the large positive correction to the masses due to in the RGB's. Also, one 
finds that the B, is typically the lightest sparticle. But most importantly, notice that one of the 
Higgs mass^, goes negative triggering electroweak symmetry breaking |25| . (The negative sign 
in the figure refers to the sign of the mass squared, even though it is the mass of the sparticles 
which is depicted.) 
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Figure 1: RG evolution of the mass parameters in the CMSSM. I thank Toby Falk for providing 
this figure. 

For given values of tan/3, Aq, and sgn{ii), the regions of the CMSSM parameter space that 
yield an acceptable relic density and satisfy the other phenomenological constraints may be 
displayed in the {nii/2,mQ) plane. In Fig. [2^, the light shaded region corresponds to that 
portion of the CMSSM plane with tan/3 = 10, Aq = 0, and /i > such that the computed relic 
density yields the WMAP value given in eq. ^TE\\ [22] . The bulk region at relatively low values of 
™i/2 a-nd TOO; tapers off as toi/2 is increased. At higher values of toq, annihilation cross sections 
are too small to maintain an acceptable relic density and Q,^h^ is too large. Although sfermion 
masses are also enhanced at large TOi/2 (due to RGE running), co-annihilation processes between 
the LSP and the next lightest sparticle (in this case the f) enhance the annihilation cross section 
and reduce the relic density. This occurs when the LSP and NLSP are nearly degenerate in 
mass. The dark shaded region has rrif < and is excluded. The effect of coannihilations is 
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to create an allowed band about 25-50 GeV wide in toq for mi/2 ^ 950 GeV, or mi/2 ^ 400 
GeV, which tracks above the m^j = contour pB] , 




mi/2 (GeV) mi/2 (GeV) 



Figure 2: The (mi/2,mo) planes for (a) tan/3 = 10 and ^ > 0, assuming = 0, = 175 GeV 

and mb(mb)g^ = 4.25 GeV. The near-vertical (red) dot-dashed lines are the contours m/i = 
114 GeV, and the near- vertical (black) dashed line is the contour m^± = 104 GeV. Also shown 
by the dot-dashed curve in the lower left is the corner excluded by the LEP bound of > 99 
GeV. The medium (dark green) shaded region is excluded by 6 — J- 57, and the light (turquoise) 
shaded area is the cosmologically preferred region. In the dark (brick red) shaded region, the 
LSP is the charged fi. The region allowed by the E821 measurement of at the 2-a level, is 
shaded (pink) and bounded by solid black lines, with dashed lines indicating the 1-a ranges. In 
(b), tan/3 = 50. 

Also shown in Fig. [2^ are the relevant phenomenological constraints. These include the 
LEP lower limits on the chargino mass: m^± > 104 GeV [57] and on the Higgs mass: m/j > 
114 GcV [55]. FeynHiggs [55] is used for the calculation of irifi- The Higgs limit imposes 
important constraints, principally on mi/2 and particularly at low tan/3. Another constraint 
is the requirement that the branching ratio for — > 57 be consistent with the experimental 
measurements |30| . These measurements agree with the Standard Model, and therefore provide 
bounds on MSSM particles [31], such as the chargino and charged Higgs bosons, in particular. 
The constraint imposed by measurements oi b sj also exclude small values of mi/2- Finally, 
there arc regions of the (toi/2, rno) plane that are favored by the BNL measurement |35] of the 
anomalous magnetic moment of the muon or — 2. Here, we assume the Standard Model 
calculation [33] of g^~2, and indicate by dashed and solid lines the contours of 1- and 2-(j level 
deviations induced by supersymmetry. 

At larger mi/2,"T-o and tan/3, the relic neutralino density may be reduced by rapid annihi- 
lation through direct-channel H,A Higgs bosons, as seen in Fig.[2Ib) [Tni[3T]. Finally, the relic 
density can again be brought down into the WMAP range at large mg (not shown in Fig. [2]), 
in the 'focus-point' region close the boundary where electroweak symmetry breaking ceases to 
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be possible and the lightest neutralino x acquires a significant higgsino component [34]. 

As seen in Fig. [21 the relic density constraint is compatible with relatively large values of 
™i/2 a-nd mo. However, all values of mi/2 and toq arc not equally viable when the available 
phenomenological and cosmological constraints are taken into account. A global likelihood 
analysis enables one to pin down the available parameter space in the CMSSM. One can avoid 
the dependence on priors by performing a pure likelihood analysis as in [3S] , or a purely x^-based 
fit as done in [311 [37] ■ Here, the results from one such analysis [35] [3S] is presented, which used 
a Markov-Chain Monte Carlo (MCMC) technique to explore efficiently the likelihood function 
in the parameter space of the CMSSM. A full list of the observables and the values assumed 
for them in this global analysis are given in [37] , as updated in [SHJ [39] . 

The best fit point is shown in Fig. [3] which also displays contours of the A^^ function in the 
CMSSM. The parameters of the best-fit CMSSM point arc mo = 60 GeV, mi/2 = 310 GcV, 
^0 = 130 GcV, tan/3 = 11 and ^ = 400 GcV, yielding the overall xVNdof = 20.6/19 (36% 
probability) and nominally Mh = 114.2 GeV [Mj. The best-fit point is in the coannihilation 
region of the (jtiq, mj^/2) plane. The C.L. contours extend to slightly larger values of mo in the 
CMSSM. However, the qualitative features of the Ax^ contours indicate a preference for small 
Too a-nd mi/2- It was found in [38| that the focus-point region was disfavored at beyond the 
95% C.L. in the CMSSM. We see in Fig. [3] that this region is disfavored at the level Ax^ ~ 8 
in the CMSSM. 




nio [GeV/c^] 



Figure 3: The Ax^ functions in the (TOo,mi/2) planes for the CMSSM. We see that the coan- 
nihilation region at low too and mi/2 is favored. 

As noted above, there are several important cosmological and phenomenological constraints 
on the supersymmetric parameter space. Improvements in sensitivity have made it possible 
for direct detection experiments [101 [S] to be competitive as well. The elastic cross section 
for X scattering on a nucleus can be decomposed into a scalar (spin- independent) and a spin- 
dependent part. Each of these can be written in terms of the cross sections for elastic scattering 
off individual nucleons. The scalar part of the cross section can be written as 

asi = — [ZU + Z)f^f , (18) 
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where rrir is the x-nuclear reduced mass and 



TTiN ^ m„ 27'^™ ^ m„' ^ ' 



q—u.d,s ^ q—c.b,t 

for — p or n. The parameters are defined by 

m^/^^' = {N\m,qq\N) = m^B^, (20) 

and the a^g contain the individual quark- neutrahno scattering cross sections, see [121 HI] 
for further details regarding the calculation of the cross section. 

The elastic scattering of neutralinos on nucleons is very sensitive to the strangeness contri- 
bution to the nucleon mass and can be characterized by the parameter, y, which is also related 
to the TT-nucleon sigma term Sttn by 

y ^ R^f R = 1 - '^n/^- ■ (21) 

where ctq is the change in the nucleon mass due to nonzero u and d masses and is estimated 
from octet baryon mass differences to be ctq = 36 MeV [45], and the latest determination 
of = 64 MeV. The effect of varying these assumptions are discussed in the context of the 
CMSSM in [331 133] . Lattice calculations are now reaching the stage where they may also provide 
useful information on Sttat [46] . and a recent analysis would suggest a lower value Sttjv Si 40 

In panel (a) of Fig.|3]the spin-independent neutralino- nucleon elastic scattering cross sections 
are shown as functions of neutralino mass for the regions of Fig. [2^ that are cosmologically viable 
(i.e., those where the upper limit on the relic density of neutralinos is respected), and are not 
excluded by constraints from colliders. Here, however, parameter values corresponding to the 
focus point at high niQ are also included. Also plotted are the limits on the spin-independent 
cross section from CDMS II gS] (solid black line) and XENONIO (solid red line), as well 
as the sensitivities projected for XENONIOO [49] (or a similar 100-kg liquid noble-gas detector 
such as LUX, dashed red line) and SuperCDMS at the Soudan Mine [SO] (dashed black line). 

There are two distinct regions in the {rriy^^cr) plane, that arising from the focus-point region 
at < 150 GeV and relatively large cr, and that from the coannihilation strip. In the 
coannihilation strip, 50 GeV < < 400 GeV, where the lower limit on is a result of 
the LEP constraint on the chargino mass, and the upper limit on corresponds to the end- 
point of the coannihilation strip for tan /3 = 10. In contrast, the end point of the focus-point 
region shown is due only to the cut-off mp < 2 TeV that has been assumed. In addition, for 
"^1/2 ^ 380 GeV in the coannihilation strip (m^ < 160 GeV), the nominal calculated mass 
of the lighter scalar MSSM Higgs boson is less than the LEP lower bound. These points are 
indicated by lighter shadings. 

The choices tan/3 = 10 and = do not yield viable direct detection cross sections that 
are completely representative of the range of possibilities within the CMSSM. Therefore, in 
Fig. |4}d, we show CMSSM spin-independent neutralino-nuclcon cross section, as obtained in 
a scan over aU CMSSM parameters with 5 < tan/3 < 55, 100 < mi/2 < 2000 GeV, GeV 
< 7710 < 2000 GeV, and — 3mi/2 < Aq < 3mi/2 [3H]- We also allow both positive and negative 
/X, except for large tan /? > 30, where convergence becomes difficult in the < case. At low 
< 300 GeV, cross sections generally exceed 10~^ pb, and the largest scalar cross sections. 
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Figure 4: The neutralino-nucleon cross sections as functions of neutralino mass for the CMSSM 
with tan/3 = 10. Also shown are upper Umits on the cross section from CDMS II [40] 
(solid black line) and XENONIO [JT] (solid pink line), as well as the expected sensitivities 
for XENONIOO (dashed pink line) and SuperCDMS at the Soudan Mine [50] (dashed black 
line) . Panel (b) shows the neutralino-nucleon cross sections as functions of neutralino mass for 
the CMSSM, with 5 < tan/3 < 55, < mi/2 < 2000 GeV, 100 GeV < mo < 2000 GeV, and 
— 3mi/2 < j4o < 3mi/2 We consider /i < only for tan/3 < 30. Taken from [48] . 



which occur for - 100 GeV, are aheady excluded by CDMS II [40] and/or XENONIO 
[41| . These exclusions occur primarily in the focus-point region at large tan/3. On the other 
hand, for rriy^ ^ 400 GeV scalar cross sections are well below 10 ~^ pb, and come from the 
coannihilation strip or the rapid-annihilation funnel that appears at large tan/3 in the CMSSM. 
The effective cross sections shown are suppressed for points with ^cdm, and there may 

be cancellations at larger that suppress the cross sections substantially. These regions of 
parameter space will not be probed by direct detection experiments in the near future [49 [ 150] . 

Finally, the frequentist analysis described above can also be used to predict the neutralino- 
nucleon elastic scattering cross section [5^. The value of a^^ shown in Fig. [S^ is calculated 
assuming a tt-N scattering a term Eat = 64 MeV. We see in Fig. [5] that values of the Xi-proton 
cross section ct^^ ^ 10~* pb are expected in the CMSSM, and that much larger values seem 
quite unlikely. The 2D function is shown in Fig. [51d. 
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